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ABSTRACT 

We investigate the nature of Lyman break galaxies (LBGs) using a chemodynamical 
model for evolution of galaxies. Our models predict an early (the first Gyr) stage of 
intense star formation in the evolution of massive spheroids which could be identified to 
the LBGs, observed at redshift ~ 3 with strong ongoing star formation. In particular, 
we are successful in reproducing the properties of the LBG DSF 2237+116 C2 with 
a model describing a young ~ L* spheroid. The comparison of the predictions of our 
models with the observations gives support to the scenario in which LBGs are the 
progenitors of present-day massive spheroids, i.e. bulges of luminous early type spirals 
or luminous elliptical galaxies. 

Key words: cosmology: observations - galaxies: elliptical - galaxies: evolution - 
galaxies: formation- galaxies: ISM - galaxies: starburst 



1 INTRODUCTION 



Colour selection techniques based on the Lyman limit break 
of the spectral energy distribution caused by neutral hydro- 
gen absorption have been used for many years in surveys for 
distant QSOs (e.g. Warren et al. 1987). Guhathakurta et al. 
(1990) and Songaila, Cowie & Lilly (1990) used this method 
to set limits on the number of star-forming galaxies at z w 3 
in faint galaxy samples. More recently, Steidel & Hamilton 
(1992, 1993) and Steidel, Hamilton & Pettini (1995), using 
this method, designed a broad band filter set (the U„G1Z 
system), which allowed them to discover a widespread pop- 
ulation of star forming galaxies at redshift z ~ 3, the Ly- 
man break galaxies (LBGs). Spectroscopic confirmation of 
their redshifts was first presented by Steidel et al. (1996), 
and WFPC2 images of select LBGs were published by Gi- 
avalisco, Steidel & Macchetto (1996). 

An important recent advance in the study of LBGs was 
the availability of the first results from a program of near- 
infrared spectroscopy aimed at studying the familiar rest- 
frame optical emission lines from H II regions of LBGs (Pet- 
tini et al 1998b, hereafter P98). The program was successful 
in detecting Balmer and [O III] emission lines in five LBGs. 
The nebular luminosities imply star formation rates (SFRs) 
larger than those deduced from the UV continuum, which 
suggests significant dust reddening. In four LBGs the ve- 
locity dispersion of the emission lines is <r em ~ 70 km s~\ 
while the fifth system has a em — 200 km s _1 . The relative 
redshifts of interstellar absorption, nebular emission, and 
Lyman a emission lines differ by several hundred km s _1 , 
a similar effect to that found in nearby HII galaxies ( Kunth 



et al 1998) indicating that large-scale outflows may be a 
common characteristic of both starbursts and LBGs. 

On the other hand, we have developed a chemodynam- 
ical model (Friaga & Terlevich 1994; Friaga & Terlevich 
1998, hereafter FT) for formation and evolution of spheroids, 
which are suspect to be the z = counterparts of LBGs 
(Steidel et al. 1996). Our chemodynamical model combines 
multi-zone chemical evolution with 1-D hydrodynamics to 
follow in detail the evolution and radial behaviour of gas and 
stars during the formation of an spheroid. The star forma- 
tion and the subsequent stellar feedback regulate episodes 
of wind, outflow, and cooling flow. The knowledge of the 
radial gas flows in the galaxy allows us to trace metallic- 
ity gradients, and, in particular, the formation of a high- 
metallicity core in ellipticals. The first ~ 1 Gyr of our model 
galaxies shows striking similarities to the LBGs: intense star 
formation, compact morphology, the presence of outflows, 
and significant metal content. We now proceed to exam- 
ine these similarities, and, in particular, to consider the im- 
plications of the recent near-infrared observations of P98. 
We demonstrate that our model supports the scenario in 
which LBGs are the progenitors of the present-day bright 
spheroids. In this paper, the SFRs, luminosities and sizes 
quoted by P98 are converted to the cosmology adopted here 
{Ho = 50 km s" 1 Mpc" 1 , q = 0.5). 
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2 LYMAN BREAK GALAXIES AS YOUNG 
SPHEROIDS 

There are several evidences in favour of the LBGs being 
the high-redshift counterparts of the present-day spheroidal 
component of luminous galaxies (Steidel et al. 1996, Gi- 
avalisco et al. 1996): their comoving space density is at least 
25 % of that of luminous (L > L*) present-day galaxies; 
the widths of the UV interstellar absorption lines in their 
spectra imply velocity dispersions of 180 — 320 km s _1 , typ- 
ical of the potential well depth of luminous spheroids; they 
have enough binding energy to remain relatively compact 
despite the very high SN rate implied by their SFRs. In ad- 
dition, the population of LBGs shows strong clustering in 
concentrations which may be the precursors of the present 
rich clusters of galaxies at a time when they were beginning 
to decouple from the Hubble flow (Steidel et al. 1998). In 
the context of Cold Dark Matter models of structure forma- 
tion, the LBGs must be associated with very large halos, of 
mass <; 10 12 Mq, in order to have developed such strong 
clustering at z ~ 3. 

Assuming a Salpeter IMF, P98 inferred from the emis- 
sion Balmer lines values for the SFR (uncorrected for dust) 
of their LBGs in the range 19-210/i5 2 M yr -1 . These val- 
ues are typically a factor of several larger than those deduced 
from the UV continuum and indicate that the correction for 
dust is typically 1-2 magnitudes at 1500 A. Dickinson (1998), 
for a large sample of LBGs, deduced from the UV continuum 
SFRs in the range 3 — 60/igQ 2 M© yr -1 . Assuming a 1 Gyr 
old continuous star formation, he used the G — 1Z colours to 
compute corrections for dust extinction to the SFR. With 
a Calzetti (1997) attenuation law, after correction for dust 
extinction, the SFR range becomes ~ 3— ~ 1500 Mq yr -1 . 
These levels of star formation are remarkably close to the 
values of the SFR exhibited in the early evolution of the 
chemodynamical models of FT. 

FT built a sequence of chemodynamical models repro- 
ducing the main properties of elliptical galaxies. The calcula- 
tions begin with a gaseous protogalaxy with initial baryonic 
mass Mg- Intense star formation during the early stages of 
the galaxy builds up the stellar body of the galaxy, and dur- 
ing the evolution of the galaxy, gas and stars exchange mass 
through star formation and stellar gas return. Owing to in- 
flow and galactic wind episodes occuring during the galaxy 
evolution, its present stellar mass is ~ 15 — 70% higher than 
Mg- Gas and stars are embedded in a dark halo of core 
radius ru and mass Mh (we set Mh = 3Mg). The models 
are characterised by Mg, r h , and a star formation prescrip- 
tion. The SFR is given by a Schmidt law vsf oc p" SF (p 
is the gas density and vsf = SFR/p is the specific SFR). 
Here we consider the standard star formation prescription 
of FT, in which the normalization of v is vq = 10 Gyr -1 
(in order to reproduce the suprasolar [Mg/Fe] ratio of gi- 
ant ellipticals), usf = 1/2, and the stars form in a Salpeter 
IMF from 0.1 to 100 M . A more detailed account of the 
models can be found in FT. Figure 1 shows the evolution of 
the SFR for models with M G in the range 5 x 10 9 - 5 x 10 11 
Mq (r h — 0.8 — 5 kpc). During the maximum of the SFR, 
the stellar velocities dispersions of these models, 55 — 220 
km s _1 , bracket the a em = 55 — 190 km s _1 range of the 
P98's LBGs. The corresponding present-day (age of 13 Gyr) 
luminosities are 0.05L* - 1AL* (-M B = 17.6 - 21.3). For 



our models, the typical range of SFR averaged over the first 
Gyr, 10 - 700 M yr -1 , reproduces well the SFRs found for 
LBGs, deduced from both the Balmer lines and the UV con- 
tinuum corrected for dust extinction. In addition, the SFR 
drops dramatically after 1.5-2 Gyr, and becomes below the 
lowest SFRs found for the LBGs. The similarity of of the 
SFRs of our models to those of LBGs allows us to identify 
the LBGs to young (<J 1 — 2 Gyr) spheroids. 

It is important to note that the moderately high SFRs 
of the LBGs seem to be difficult to conciliate with the 
predictions of the simplistic one-zone (or monolithic) mod- 
els of formation of elliptical galaxies for supra-L* systems. 
The monolithic models of formation of early-type galaxies 
have been worked out in the early 1970's (e.g. Larson 1975) 
and are succesful at reproducing the supra-solar [Mg/Fe] of 
bright ellipticals (Matteuccci & Tornambe 1987; Hamann & 
Ferland 1993), but the required short star formation time 
scale (~ 10 yr) implies extremely high SFRs during the 
formation of L > L* ellipticals. As a matter of fact, in 
the one-zone model, a gaseous protogalaxy with 5 x 10 10 
M , would have a peak SFR of ~ 5000 M yr _1 , and a 
present-day Mb = —21.1. At least at redshift 3 ^ z ^ 3.5, 
such SFR is excluded by the properties of the population 
of LBGs. By contrast, in the chemodynamical model, the 
metallicity and abundance ratios of the central region of the 
young elliptical are explained with no need for all the galaxy 
having a global starburst coordinated with the central star- 
burst, which avoids the excessively high SFRs of the one- 
zone model. The most massive model here (Mg = 5 x 10 11 
M ; present-day M B = -21.3) has a peak SFR of 1050 
M yr _1 , consistent, after correction for dust extinction, 
with the highest SFRs derived from the UV continuum of 
LBGs (Dickinson 1998). Note that, as we show below, be- 
cause the observed rest-frame UV colours limit the amount 
of dust extinction to ~ 3 mag at most, we cannot evoke dust 
to hide a 5000 Mq yr" 1 starburst as a LBG at z ~ 3. 

HST optical imaging, which probes the rest frame UV 
between 1400 and 1900 A, has revealed that the LBGs 
are generally compact, with a typical half-light radius of 
1.4-2. Ih^o kpc (Giavalisco et al. 1996). The observed LBGs 
do not seem to have disk morphology, with the exception of a 
few objects without central concentration. In addition, some 
objects have a light profile following a r 1 ^ 4 law over a large 
radial range, which supports the identification of LBGs to 
young spheroids. Near infrared imaging have yielded half- 
light radii in the range 1.7 — 2.3/iJg 1 kpc (P98). The sim- 
ilarity of the near-infrared sizes to those obtained by the 
HST suggests that the optical morphology follows the UV 
morphology. As shown in the next section, the compact ap- 
pearance of the LBGs, both in the UV and in the optical, is 
reproduced by our young spheroid models. 

Note that, due to the strong fading of surface brightness 
with redshift (oc (1 + z)~ 4 ), the outer parts (r <; 10 kpc) of 
the galaxy with milder star formation rates (vsf ~ 1 Gyr -1 
or less) would be missed in high redshift observations. The 
difficulty in observing the outer regions of the galaxy would 
only be compounded if there is some dust extinction. There 
is an analogy between the LBGs and nearby HII galaxies, in 
which we are observing only the brightest part of the galaxy, 
superposed on much more extended low surface brightness 
object, when deeper expositions are made available (Telles 
& Terlevich 1997; Telles, Melnick & Terlevich 1997). Addi- 
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tional support to the LBG-starburst connection comes from 
the fact that the LBGs in the P98 sample fall on the ex- 
trapolation to higher luminosities of the correlation I/ H/ 3 — o 
found for local H II galaxies by Melnick, Terlevich, & Moles 
(1988) (Terlevich 1998). 



3 DSF 2237+116 C2, A YOUNG ~ L* 
SPHEROID? 

It is of interest to compare the predictions for our models 
with the observational data of DSF 2237+116 C2, the most 
massive LBG (the LBG with the largest a em ) in the P92 
sample. The properties of this object are successfully de- 
scribed by the fiducial model of FT (M G = 2 x 10 11 M 
and Th = 3.5 kpc ). Its present-day stellar mass, 2.4 x 10 11 
M , corresponds to Lb = 0.7L* , which allows us to iden- 
tify DSF 2237+116 C2 to an ~ L* spheroid seen during its 
early evolution, characterised by intense star formation. For 
the fiducial model, Figure 1 shows the evolution of the SFR 
within several radii. The initial stage of violent star forma- 
tion lasts ~ 1 Gyr, and exhibits a maximum SFR of ~ 500 
Mq yr" at 0.6 Gyr. After the galactic wind is established 
(at t = 1.17 Gyr), the SFR plummets and practically all star 
formation within 10 kpc is concentrated inside the inner kpc. 
The late central star formation, characterised by a moderate 
SFR (~ few M^yr" 1 ), is fed by a cooling flow towards the 
galactic centre. The stagnation point separating the wind 
and the inflow moves inwards until it reaches the galactic 
core at t = 1.8 Gyr, when a total wind is present through- 
out the galaxy. After this time, indicating the end of the 
star-forming stage, only very small levels of star formation 
are present in the galaxy. The early stage of star formation 
during which the stellar body of the galaxy is formed (the 
stellar mass reaches 50% of its present value at t = 3.9 x 10 8 
yr), resembles the LBGs. The average SFR during the first 
Gyr, 328 M yr" 1 , is very similar to the SFR of 210h^ M Q 
yr" 1 of DSF 2237+116 C2 inferred from its H/3 luminosity. 
In addition, the SFR is concentrated in the inner 2-3 kpc, 
which gives to our model galaxy the compact appearance 
typical of LBGs. 

Figure 1 also shows Lisoo, the luminosity at 1500 A, 
which allows a more direct comparison with the imaging 
data. Note that our models reproduce the compact appear- 
ance of LBG, the light being concentrated in the inner ~ 3 
kpc until the maximum of the SFR and in the inner ~ 2 
kpc after that time. The luminosities predicted during the 
first Gyr are around 3 x 10 42 erg s" 1 A -1 . This value is 
higher than the observed L1500 of4.1xl0 41 /i5 1 ergs" 1 A" 1 
found for DSF 2237+116 C2. Note that the oc (1 + z)" 4 
dimming of the surface brightness with the redshift makes 
it difficult to detect the outer regions of the galaxy. How- 
ever, considering the UV emission inside a projected ra- 
dius of 10 kpc, reduces only slightly the UV luminosity 
(L 1B0O (r < lOkpc) = 2.5 x 10 42 erg s" 1 A" 1 ). On the other 
hand, a simple comparison between the SFR deduced from 
the H/3 line, assuming that the extinction at the H/3 is neg- 
ligible, and the SFR deduced from the UV continuum, in- 
dicates for DSF 2237+116 C2 a correction factor for dust 
between 7 and 48 (P98). These very high correction fac- 
tors should not be taken at face value, since this simplistic 
approach furnishes some unphysical results, such as nega- 
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Figure 1. Top panel: evolution of the SFR for the fiducial over 
the whole galaxy and inside several radii. Also given the predicted 
SFRs averaged over the time spans 0-0.5, 0.5-1, 1-1.5, 1.5-2, 2-2.5, 
and 2.5-3 Gyr for several models (symbols labelled by the model 
Mq). Middle panel: evolution of the 1500 A luminosity over the 
whole galaxy and inside several projected radii. The luminosi- 
ties have been calculated from the models of Bruzual & Chariot 
(1998), for a Salpeter IMF from 0.1 to 100 Mq and metallici- 
ties from Z = 0.0001 to Z = 0.1. Lower panel: evolution of the 
unreddened G — TZ colour for the fiducial model, over the whole 
galaxy and inside a projected radius of 10 kpc (dotted and dot- 
dashed lines, respectively). Also shown (with the same line styles) 
the amount of reddening A1500 needed for G — 1Z of the model 
to match the observed G — TZ. In all three panels the dashed 
horizontal lines represent the respective observed quantities for 
DSF 2237+116 C2: SFR deduced from the H/3 emission; L1500 
uncorrected for dust; G — 1Z colour. In the middle panel, the two 
thin horizontal lines denote the observed L1500 with a correction 
for dust extinction deduced from the G — IZ colour as in P98 
(see text), and the dotted line indicates the observed L1500 cor- 
rected for dust extinction using the A1500 value over the whole 
galaxy shown in the lower panel. In the upper panel, the thin long 
dashes on the right denote the SFRs deduced from H/3 emission 
for the <r 70 km s _1 LBGs of P98; the highest value corre- 
sponds in fact to an upper limit for Q0000-263 D6, since for this 
object (at z em = 2.966), H/3 is outside the K-band, and the SFR 
was deduced from the [O IIL1A5007 luminosity, assumimg H/3/[0 
IILU5007 < 0.5. 
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tive extinctions for some objects. It would be interesting to 
consider a dust extinction index based on the UV part of 
the spectrum, the most easily accessible to observations of 
LBGs. The effect of dust is to flatten the spectrum, and the 
colour G — TZ provides a reliable measure of the UV slope 
(at z ~ 3, the effective redshifts of the two filters, 4740 
and 6850 A, respectively, are translated to 1190 and 1710 
A). The comparison of the observed (G — TZ) bs colours to 
the (G — lZ) ca i c colours predicted by an unreddened contin- 
uous star formation model with absorption by the Lyman 
a forest, allowed P98 to deduce dust correction factors be- 
tween ~ 1 and ~ 10 for the UV luminosities of the LBGs in 
their sample. In the case of DSF 2237+116 C2, a value of 
iisoo = 3.9(1.8) x 10 42 erg s -1 A -1 is obtained after a cor- 
rection for dust extinction assuming a Calzetti attenuation 
law and a continuous 10 7 (10 9 ) years old star formation. 

In view of the importance of the G — TZ colour in check- 
ing for star formation and estimating the dust extinction, 
Figure 1 also shows (G-TZ) ca ic predicted for DSF 2237+116 
C2 (z = 3.317) by the fiducial model, obtained as follows: 
in the first place, the integrated SED is calculated for sev- 
eral apertures, using the Bruzual & Chariot (1998) mod- 
els; then the SED is redshifted to z — 3.317, reddened 
by the Lyman a forest opacity (Madau 1995), and con- 
volved with the filter transmission curves. Finally, when 
(G — TZ)caic is bluer than the G — TZ colour of the galaxy 
((G — 7?.) 6 3 =1.13), we calculate, assuming a Calzetti atten- 
uation curve, the value of A1500 needed to match (G — TZ) bs- 
Since G — 1Z becomes redder with time, we can use the condi- 
tion (G — TZ) ca ic<(G — TV) bs to set an upper limit in the age 
of the galaxy, beyond which ^1500 becomes formally nega- 
tive. This limit is 1.00 Gyr, for an aperture r < 10 kpc, and 
1.52 Gyr, if the aperture encompasses the whole galaxy. The 
predicted colours are bluer for the larger aperture because: 
1) metallicities are typically ~ 0.1 solar for r > 10 kpc, im- 
plying bluer colours for the star population; and 2) there is 
some star formation in the outer parts of the galaxy as the 
gas driven by the galactic wind is compressed on its way out 
of the galaxy. At the peak of the SFR, A1500 reaches w 2.15, 
within the range A1500 = 1-58 — 2.44 deduced by P98 for a 
continuous star formation lasting from 10 9 to 10 7 yr. Figure 
1 also shows the observed value of L1500 corrected for dust 
extinction using the time- dependent value of A1500 obtained 
as above, and also the values corrected as in P98. The agree- 
ment with the predictions of our models both for the galaxy 
as a whole as for the inner 10 kpc is excellent. Therefore, if 
our model galaxy were at a redshift ~ 3, it would be easily 
seen as an LBG. 

In order to explore the recent availability of infrared 
imaging, tracing the rest-frame optical light, Figure 2 shows 
the blue luminosity of the fiducial model inside several pro- 
jected radii. The similarity of rest-frame optical and UV 
sizes, indicated by the optical and near-infrared observa- 
tions, is reproduced by the predictions of our model: the 
half-light radii at the maximum of SFR, at 1500 A and in the 
blue band, are 1.64 and 1.51 kpc, respectively. It is useful, 
due to the possibility of missing light from the outer parts of 
the galaxy, to consider the half-light radii with respect only 
to the inner 10 kpc of the galaxy. In this case, the half-light 
radii at the SFR peak are 1.46 and 1.44 kpc, for 1500 A and 
blue light, respectively. Therefore, the optical morphology 
follows the UV morphology, and the galaxy remains com- 



pact in the optical band. Note, however, that the light does 
not trace the mass. At the maximum of SFR, the half-mass 
radius (=7.5 kpc) is much larger than the half-light radius. 
The star formation does not follow the stellar mass, but in- 
stead it is regulated by the gas flows (e.g., the star formation 
within the inner kpc is fed by the cooling flow towards the 
galaxy centre). The star formation is not coordinated along 
the galaxy: vsf in the inner kpc reaches several xlO Gyr -1 , 
whereas the vsf averaged over the whole galaxy is slightly 
larger than 1 Gyr - . In view of this, estimating the mass 
from the half-light radius will seriously underestimate the 
galaxy mass. P98 were suspicious of having underestimat- 
ing the mass of DSF 2237+116 C2 (the value they derive 
is 5.5 x 10 10 Mq). Here we quantify their suspicion, sug- 
gesting that the mass underestimate could be a factor 4-5. 
In fact, at the SFR peak, our model predicts not only half- 
light radii (whatever their definition) that are very similar 
to the 1.7/igo 1 found for DSF 2237+116 C2, but also a stel- 
lar velocity dispersion of 179 km s _1 , essentially identical 
to the observed cr em = 190 ± 25, whereas the stellar mass of 
our galaxy model is 2 x 10 11 M Q at this time. 

The metal lines in the spectra of LBGs, with origin 
in stellar photospheres, interstellar absorption, and nebular 
emission, indicate metallicities anywhere between 0.01 solar 
and solar (Steidel et al. 1996). On the other hand, the strong 
correlation between the UV spectral index and metallicity in 
local starbursts would suggest a broad range in metallicity 
from substantially subsolar to solar or higher (Heckman et 
al. 1998) . In order to make predictions on the metal content 
of LBGs, Figure 2 also shows the average metallicity of the 
stellar population, inside several spherical zones. The inner 
region reaches solar metallicities (at 1.11 x 10 8 and 1.56 x 
10 s for the inner kpc and for the 1 < r < 2 kpc region, 
respectively) much earlier than the maximum in the SFR. 
Therefore, when the galaxy becomes visible as a LBG (i.e. as 
a star- forming galaxy), its metallicity inside a typical half- 
light radius (~ 1.5 kpc) will be solar or suprasolar. On the 
other hand, substantial abundance gradients are built up. 
The metallicity approaches 3 Zq in the inner kpc, while it 
is typically ~ 0.1 Z Q for r > 10 kpc. 

Other important success of our models is the predic- 
tion of important outflows during the stage of intense star 
formation, which could account for the outflow at a veloc- 
ity of 500 — 1000 km s -1 in the interstellar medium of DSF 
2237+116 C2 suggested by the relative velocities of the Ly- 
man a emission lines and of the interstellar absorption lines. 
As a matter of fact, following the maximum of the SFR, an 
outflow appears at the intermediate radii, between 2 and 10 
kpc. As we can see from Figure 2, once the outflow in the 
intermediate region is established, outflow flows velocities of 
500 - 1000 km s -1 are achieved for t ~ 1 Gyr. After 1.17 
Gyr, when the outflow reaches the galaxy tidal radius (i.e. 
the onset of the galactic wind), the wind velocity increases 
up to about 1900 km s _1 . However, during the late galactic 
wind stage, the density in the outflowing gas drops dramat- 
ically, making it difficult to obtain any signature of the out- 
flow via interstellar absorption lines and emission lines. The 
flow structure is complex, because at inner radii there is a 
highly subsonic (inflow velocity 10 km s -1 ) cooling flow, 
and through the outer tidal there is infall of low density gas 
proceeding at 60 km s~ . Therefore, the high density, high 
velocity outflowing gas in the intermediate region just after 



© 1998 RAS, MNRAS 000, 000-000 



Lyman break galaxies and spheroids 5 




100 1000 10 4 10 5 

r(pc) 

Figure 2. Top panel: evolution of the blue luminosity of the fidu- 
cial model inside several projected radii. Middle panel: evolution 
of the average mctallicity of the stellar population, within several 
spherical zones, illustrating the presence of metallicity gradients 
and the time scales for chemical enrichment. The solar abun- 
dances are taken from Grevesse & Anders (1989). Lower panel: 
evolution of the velocity profile, showing the onset of the outflow 
at intermediate radii. 

the peak in the SFR explains the large scale outflows with 
velocities of ~ 500 km s~\ deduced from the relative red- 
shifts of the interstellar absorption and Lyman a emission 
lines, which are a common feature of LBGs (P98) . 

The success of our model for a ~ L* spheroid or ellip- 
tical galaxy in reproducing several properties of the LBG 
DSF 2237+116 C2 gives additional support to the scenario 
in which LBGs are the progenitors of present-day bright 
spheroids. High angular resolution spectroscopy will in the 
future provide important information regarding the velocity 
field and angular momentum of LBGs and help us to discern 
if they are young bulges or young ellipticals. 



4 DISCUSSION 

The agreement of the fiducial model with the properties of 
DSF2237-C2 suggests that the mass range of LBGs does in- 
clude present day ~ L* objects. Note that the present model 



not only accounts for this particularly massive LBG but also 
successfully predicts the properties of the ensemble of the 
LBGs, within the scenario in which they are the progenitors 
of the present day spheroids with 0.1L* Lb Si L* . Our 
models also reproduce the main properties of the four LBGs 
with lower <7 em 's in P98. This is illustrated in Figure 1, in 
which the SFRs deduced for these LBGs are similar to those 
of models with M G = 10 10 - 5 x 10 10 M , for t <5 1 - 1.5 
Gyr (present day -M B = 18.4 - 19.7). 

These models chosen because they exhibit during the 
period 0.2 < t < 1.5 Gyr (the lower limit on time garan- 
tees that a significant stellar component has already been 
formed, and for times later than the upper limit, THE SFR 
has probably decreased below levels typical of LBGs) the 
stellar velocity dispersion coincides with the values of a cm 
of the 4 low (T em LBGs of P98 (which are in the range 
55 ±15 -85 ±15 kms^ 1 ). 

One of the central aspects of our modelling is that it 
follows in detail the impact of gas flows on the early evolu- 
tion of the galaxies. Besides the importance of galactic winds 
in galaxy evolution, as already highlighted in the pioneer- 
ing work of Larson (1974), cooling flows also play a central 
role in galaxy evolution — feeding a central AGN hosted 
in the galatic core, building up metallicity gradients (FT), 
and maintaining a moderate level of star formation in the 
inner regions of the galaxy at late times, i.e. when the major 
stellar population of the elliptical galaxy has already been 
formed (Jimenez et al. 1998). 

As a matter of fact, the flow structure is complex, ex- 
hibiting, for instance, during a considerable span of the 
galaxy evolution a partial wind, with inflow in the inner 
parts of the galaxy and outflow in the outskirts of the galaxy. 

Moreover, the flow structure varies with time, and the 
same star-forming galaxy can exhibit a variety of flow pro- 
files, depending on the evolutionary stage being picked up 
by the observation. As can be seen from Figure 2, in the 
fiducial model the outflow does not occur during the whole 
period of intense star formation. Outflow velocities of ~ 500 
km s _1 are achieved only ~ 0.3 Gyr after the maximum of 
the SFR and of ~ 1000 km s" 1 ~ 0.4 Gyr after the maxi- 
mum. The delay between the maximum of the SFR and the 
onset of the outflow reflects the time needed for the energy 
input by SNe into the ISM to overcome the gravitational 
binding energy of the gas. It is possible the observation of 
LBGs, i.e. with a high SFRs, both in the phase of outflow 
and before the onset of the outflow. For earlier times, there 
are global inflows, reaching velocities of up to a few 100 
km s _1 . Therefore, we expect a large dispersion in the rela- 
tive redshifts of the interstellar absorption, nebular emission 
and Lyman a emission lines of LBGs. As we discuss below, 
this seems to be the case. 

We can see from Figure 3 that for the models with Ma 
in the range 10 1() - 5 x 10 10 M Q , which describe well the 
four a em ~ 70 km s _1 in P98, the evolution of radial flows is 
qualitatively similar to that of the fiducial model. The main 
difference is that the outflow happens earlier, and, once the 
outflow is established, velocities higher than 1000 km s _1 
are reached faster. This is a result of the shallower potential 
well of these galaxies. Note however, that the final wind 
velocities are somewhat lower than in the fiducial model. 

Assuming that the Balmer and [O III] emission lines 
are at the galaxy systemic redshift, the velocity shifts of 
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Figure 3. Evolution of the gas velocity profiles of low mass 
models aimed at describing the a as 70 km s _1 galaxies of P98. 



the interstellar absorption, nebular emission, and Lyman a 
emission lines found by P98 indicate that large velocity fields 
arc a common feature of LBGs. It is important to note that 
exactly the same result has been found in nearby HII galax- 
ies (Kunth etal 1998). In all cases where Lyman a emission 
was detected, the line peak is shifted by « 1000 km s^ 1 rel- 
ative to the metal absorption lines. In two out of three cases 
(Q0000-263 D6 and B2 0902+343 C6) the nebular lines are 
at intermediate velocities, with Q0000— 263 D6 exhibiting a 
conspicuous P-Cygni profile. The most ready interpretation 
of these characteristics is the presence of large scale outflows 
with velocities of ~ 500 km s _1 in the interstellar media of 
the galaxies observed: the Lyman a emission is suppressed 
by resonant scattering and the only Lyman a photons es- 
caping unabsorbed in our direction are those back-scattered 
from the far side of the expanding nebula, whereas in absorp- 
tion against the stellar continuum we see the approaching 
part of the outflow. Within this scenario, the relative ve- 
locities of the three line sets of these LBGs are consistent 
with our model, since velocities of ~ 500 km s _1 are easily 
reached in the low mass models, at a galaxy age of 0.8 — 0.9 
Gyr. 

However, this simple symmetric picture probably does 
not account for all the variety of possible situations, since 



in the third LBG with Lyman a emission (the high a em 
DSF 2237+116 C2) H/3 and [O III] emission are apparently 
at roughly the same velocity as the absorption lines, even 
though Lyman a emission is redshifted by w 1000 km s _1 . 
Also this case could be explained by our models with the 
Lyman a emission originating in the high velocity shell re- 
ceding in the back of the galaxy (in the fiducial model, at 
t — 1.1 Gyr this shell is at r « 10 kpc with a velocity of 1250 
km s _1 ), while the interstellar absorption lines could arise in 
the few central kpc, with much lower expansion velocities. 
See the interesting discussion regarding escape of Lyman a 
photons in HII galaxies by Kunth et al (1999). 

In the P98 sample, 2 in 5 objects do not show Lyman a 
emission. In one of these systems (Q0201+113 B13) the in- 
terstellar absorption lines are redshifted by 250 km s with 
respect to the Ha emission (this object is at z ~ 2.2 and, as a 
consequence, Ha is observed instead of H/3). If this velocity 
difference is real (no error is quoted by P98 for the rela- 
tive velocity of this object), within the scenario depicted by 
our model, we could be observing the galaxy during its early 
global inflow. The remaining object in P98 sample, also with 
no Lyman a emission, Q0201+113 C6, shows a 3200 km s _1 
difference between emission and absorption lines, too large 
to be accounted for any for our models. As a matter of fact, 
this difference is so large that the line set identified as inter- 
stellar absorption could be in reality an intervening absorp- 
tion line system. We would like to point that our predictions 
regarding the importance of the gas velocity field in the es- 
cape of Lyman a photons are similar to those of Kunth and 
collaborators for nearby HII galaxies (Kunth et al. 1998, 
1999). 

It is of interest to compare the results of the present 
model with the predictions for LBGs of the semi-analytic 
models of Baugh et al. (1998) and Mo, Mao, and White 
(1998) which are based on disk formation models. These 
models are simpler than the present model and, therefore, 
their predictions are more limited. For instance, Mo et al. 
(1998) do not discuss gas flows. However, in principle, their 
model could allow for some global characterization of the 
infall associated with the formation of a disk. Note that our 
prediction of outflows allows one to distinguish the present 
model from disk models, since disk models include infall but 
do not exhibit outflows. 

One success of the models of Baugh et al. (1998) and Mo 
et al. (1998) is the good description of the clustering proper- 
ties of the LBGs within the framework of the most popular 
hierarchical models for structure formation. Note these suc- 
cess refers in fact to the clustering of halos, independently 
if they host a spheroidal or a disklike star-forming central 
galaxy. Mo et al. (1998) claims to explain also the moder- 
ate SFRs, the small sizes and the velocity dispersions of the 
LBGs. However, as pointed out by Mo et al. themselves, 
their calculation of the SFR is very sketchy. In addition, 
their correct prediction of the compact size of the LBGs is a 
consequence of identifying the LBGs with low angular mo- 
mentum objects, and so with small size for their mass. In 
connection to this last point, one of the parameters of the 
individual galaxies in Mo et al. (1998) is A, the spin param- 
eter of the halo, and for systems with A smaller than some 
critical A C rit, the gas cannot settle into a centrifugally sup- 
poted disk without first becoming self-graviting. The final 
configuration is probably spheroidal rather than disk-like. 
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In view of this, Mo et al. (1998) admits that a sizeable frac- 
tion of LBGs could be spheiroids (see their Section 3.7), 
although their model was initially designed for the forma- 
tion of disks. It is possible the existence of a population of 
disk objects among the LBGs, although the morphological 
studies of LBGs suggest that disks are a minority in the 
LBG population. HST imaging (filters F606W to F814W, 
probing the rest frame UV range 1400 - 1900 A) has shown 
that the LBGs do not seem to have disk morphology, with 
the exception of a few objects without central concentration, 
for which an exponential profile provides a good fit to their 
surface brightness distribution (Giavalisco et al. 1996). In 
addition, some objects have a light profile following a r 1 ' 4 
law over a large radial range. Near-infrared HST NICMOS 
imaging (sensitive to the rest-frame optical light) provides 
a similar picture (Giavalisco, private communication). Fu- 
ture imaging observations would be very useful to establish 
which fraction of the LBGs are disk-like systems. 

With respect to the velocity dispersions of the LBGs, 
Mo et al. (1998) reproduces the median values of the LBGS, 
but runs into some problems with the a ~ 200 km s _1 
of DSF 2237+116 C2. since their predicted stellar velocity 
dispersions are typically around 70 km s -1 for their Qq — 1 
cosmology and 120 km s" 1 for their fio = 0.3 flat cosmology. 
In the distribution of velocity dispersions predicted by Mo et 
al. (their Figure 8), the highest velocity dispersion is ~ 170 
km s _1 for the f2o = 1 cosmology and, even for the f2o = 0.3 
flat cosmology (which predicts higher velocity dispersions), 
values of a ~ 200 km s" 1 are very rare. 

One interesting possibility offered by our modelling of 
the G — 71 colors is obtaining constraints on the age of the 
LBGs within the scenario of the young spheroid model. One 
minimal constaint is that the galaxy could not be older than 
~ 1.5 Gyr, otherwise G — 1Z would be redder than the ob- 
served. If, in addition to that, we require that the SFRs 
predicted for our model be consistent with the SFR derived 
from the H/3 emission, SFR-^p, the age of the LBG is con- 
strained to be not older than ~ 1 Gyr. In this case, the 
predicted values for ylisoo, 2.9 — 1.5 (for t from 0.1 to 1 
Gyr), are comparable to those obtained by P98, because the 
evolution of the SFR of the fiducial model, with the con- 
straint t ;S 1 Gyr and the strong decrease of the SFR be- 
yond this time, includes their simple models of continuous 
star formation lasting for 10 7 yrs and 10 9 yr. On the other 
hand, P98 estimated a value of ^lisoo = 4.2, from comparing 
the SFRs deduced from H/3 and UV luminosities, assuming 
a Calzetti attenuation law and continuous star formation 
for 10 9 yrs. This value would represent a discrepancy with 
the values of A1500 predicted by our models, were it not 
for the fact that the simple comparison of H/3 and UV lu- 
minosities could lead to unreliable estimates of A\soo, as 
illustrated by some unphysical results, e.g. negative extinc- 
tions found by P98 for some objects. In determining Aisoo 
from the SFR^p/ SFRuv ratio, one should be aware of the 
uncertainties in converting L1500 to SFR and that the long 
wavelength baseline involved in this estimation increases the 
uncertainties of assuming one particular extinction law. 

Finally, the relatively high metal abundances obtained 
by our models, could be expected with some hindsight, if 
we had considered the high star formation rates derived for 
the LBGs, together with the continued (for periods ~ 1 Gyr 



long) star formation favoured in previous works (e.g. Stei- 
del et al. 1996) . One of the important consequences of these 
relativelly high metallicities coupled to the large outflow ve- 
locities is that the halo of the galaxy will be enriched in 
metals in ~ 10 8 yr, and the circungalactic environment out 
to distances of several hundreds of kpc of the galaxy will 
be contaminated in metals in ~ 1 Gyr or less. This fast 
chemical enrichment mechanism for the galactic halo and in 
intergalactic medium could explain the chemical abundances 
of quasar absorption line systems, taken as a probe of the 
gaseous galactic halo and of the intergalactic medium as for 
instance, in Lyman limit systems (Viegas & Friaga 1995) and 
in the Lyman a forest (Friaga, Viegas & Gruenwald 1998). 
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